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Abstract 

Small  sodium  clusters  consisting  of  1  to  40  atoms  are  described  as 
spheres  of  interacting  homogeneous  electron  gas  (jellium  model).  The 
static  electronic  polarizability  is  calculated  using  self-consistent 
density-functional  methods.  An  excellent  agreement  with  recent 
experimental  results  is  obtained. 


PACS  (lumbers:  71.45.Gm,  36.40^ 


Recently  Knight  et  al 1  have  reported  on  measurements  of  static 
electronic  polarizability  of  small  sodium  clusters  consisting  up  to  40 
atoms.  The  size  dependence  of  the  polarizability  is  connected  with  the 
electronic  structure  of  the  cluster  and  shows  similar  (shell)  structure  as 
found  in  earlier  measurements  of  abundance  spectra2  and  ionization 
potenti al . 3 

Theoretical  ab  initio-type  calculations  of  the  ground  state  electronic 
and  atomic  structure  of  small  metal  clusters  have  been  restricted  to 
clusters  consisting  usually  less  than  10  atoms"* ,  and  the  polarizability  has 
been  calculated  only  for  diatomic  molecules.5  However,  several  model 
calculations  assuming  a  spherical  symmetry  have  been  presented  for  small 
metal  clusters  (see  ref.  6  and  references  therein).  Typically,  the  metal 
sphere  is  described  by  non-interacting  or  interacting  electrons  moving  in  a 
spherically  symmetric  effective  potential  in  a  way  analogous  to  the  nuclear 
shell  model.  In  a  recent  paper'0  we  have  calculated  static  and  dynamical 
polarizabilities  for  small  jellium  spheres  corresponding  to  the  density  of 
metallic  lithium.  The  purpose  of  this  communication  is  to  do  a  more 
complete  study  of  the  static  polarizabi 1 i ty  of  jellium  clusters  corres¬ 
ponding  the  density  of  sodium  and  compare  the  results  with  very  recent 
experiments  of  Knight  et  al.1  A  similar  calculation  for  closed  shell 
sodium  clusters  have  been  reported  earlier  by  Ekard.  In  the  case  of 
sodium  the  comparison  of  jellium  results  to  the  experimental  ones  is 
meaningful  and  interesting  because  (i)  the  jellium  model  gives  good  results 
for  the  surface  energy8  and  vacancy  formation  energy’  of  sodium  metal 
(meaning  that  the  pseudopotential  corrections  are  small9)  and  because  (ii) 
the  stability  of  sodium  clusters  (magic  numbers)  have  been  successfully 
explained  in  terms  of  the  jellium  model. 


In  the  jellium  model  the  metal  cluster  is  described  in  terms  of  inter¬ 
acting  electron  gas  which  moves  in  the  external  electrostatic  potential 
provided  by  a  spherical  homogeneous  positive  background  charge  distribution 
n+(r)  =  no9(R-r),  where  n0  is  the  average  conduction  electron  density  of 
the  metal  and  R  is  the  radius  of  the  cluster  (R  =  N1/3RWS,  N  being  the 
number  of  atoms  in  the  cluster  and  Rws  the  bulk  VJigner-Seitz  radius). 

The  electronic  structure  and  the  (dipolar)  response  function  needed  for  the 
polarizability  determination  are  calculated  using  self-consistent  density 
functional  techniques.  For  details  of  the  method  we  refer  to  an  earlier 
publ ication.6  The  formulation  we  use  in  calculating  the  response  function 
is  strictly  valid  only  for  clusters  with  closed  electron  shells.  For  open 
shell  clusters  we  use  spin-independent  formalism  and  energy  minimizing 
fractional  occupation  numbers  for  the  degenerate  energy  levels.  This 
method  is  found  to  give  good  results  for  free  atoms. b 

The  results  are  compared  to  the  experimental  results1  in  Fig.  1.  The 
polarizability  of  the  cluster  is  normalized  dividing  by  the  polari¬ 
zability  of  a  free  atom  and  by  the  number  of  atoms  in  the  cluster.  The 
theoretical  value  (131  Oq  )  have  been  used  for  the  atomic  pol ari zabi 1 i ty 

ctj .  This  is  actually  very  close  to  the  polarizability  of  1  electron  jelli¬ 
um  cluster  (136  a^)  and  the  use  of  the  latter  would  have  given  equally 

good  overall  agreement  between  the  theory  and  experiment.  The  opening  of  a 
new  shell  causes  an  increase  in  the  polarizability  of  N  =  3,  9,  19,  21,  and 
35.  The  same  structure  is  seen  also  in  the  experimental  values.  The 
agreement  between  theory  and  experiment  is  at  best  for  closed  shell  clus¬ 
ters  il  =  2,  8,  18,  20,  34,  and  40.  For  open  shell  clusters  the  experimen¬ 
tal  values  are  slightly  above  the  theoretical  ones. 
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The  excellent  agreement  displayed  in  Fig.  1  indicates  that  many  prop¬ 
erties  of  the  clusters  can  indeed  be  understood  in  terms  of  electron  levels 
in  a  spherical  potential.  This  is  rather  surprising  since  ab  initio  calcu- 
lations  have  shown  that  the  actual  structures  of  the  smallest  clusters 
(N  <  6)  are  planar,4*10  i.e.  seemingly  far  from  spherical.  However,  the 
planar  structure  is  the  result  of  the  relaxation  of  the  atomic  position 
according  to  the  electronic  distribution  of  non-spherical  valence  electron 
orbitals  (e.g.  p-orbitals)  as  explained  in  detail  by  Martin  et  al4,  and 
does  not  necessarily  mean  that  the  valence  electron  wave  functions  would  be 
very  different  from  those  of  a  spherical  potential.  The  delicate  relation¬ 
ship  between  the  orbital  occupancy  and  structure  will  induce  more  structure 
to  the  a,j/Nai  curve  (for  small  M)  than  seen  in  the  jell  i  urn  model.  For 
example,  the  weak  minimum  at  N  =  4  in  the  polarizability  and  the  corres¬ 
ponding  maximum  in  the  ionization  potential4  is  related  to  the  structure  of 
the  cluster  and  cannot  be  predicted  by  the  jel 1 i urn  model . 

The  absolute  value  of  the  polarizabi 1 i ties  from  our  calculations  are 
about  20 j  below  the  experimental  ones,  i.e.  the  calculated  «i  is  smaller 
than  the  experimental  one.  This  atomic  pol arizabi 1 i ty  was  calculated  using 
exactly  the  same  density  functional  method  as  for  the  pol arizabi 1 i ties  of 
the  jellium  clusters.  The  fact  that  the  calculated  atomic  polarizability 
is  within  4't  of  the  pol  ari  zabi  1  i  ty  of  one-electron  jellium  cluster  indi¬ 
cates  that  the  main  discrepancy  in  the  absolute  values  of  the  pol ari zabi 1 i - 
ties  of  small  clusters  may  not  be  the  use  of  the  jellium  model  but  the 
approxi tna ti ons  made  in  solving  the  mariy-body  problem.  The  use  of  spin- 
dependent  formalism  would  change  the  absolute  values  of  open  shell  clusters 
and  produce  additional  kinks  in  the  jellium  curve  at  h-values  corresponding 
to  half-filled  shells  (due  to  the  spin  population  determined  by  Hund's 


rules).  In  real  metal  clusters  the  magnetism  would  be  connected  with  the 
actual  atomic  arrangement  and  reduced  from  the  jellium  value.10  The  large 
magnetic  moments  of  jellium  clusters  might  then  lead  to  overestimation  spin 
effects.  The  jellium  model  can  be  refined  by  taking  into  account  the 
difference  in  the  spherically  averaged  pseudopotential  between  a  real  metal 
and  jellium  or  by  using  the  (related)  pseudojellium  model.4  These 
refinements  would  not  change  the  qualitative  size  dependence  of  the 
polarizability  and  in  the  case  of  sodium  also  the  change  in  the  absolute 
values  would  be  very  small. 

In  conclusion  we  have  calculated  the  static  electronic  polarizability 
for  sodium  clusters  using  the  jellium  model  and  density  functional  tech¬ 
niques.  The  results  are  in  good  agreement  with  the  experimental  results 
giving  further  support  to  the  ideas  tnat  many  electronic  and  structural 
properties  of  sodium  clusters  can  be  understood  in  terms  of  si ngl e-parti  cl s 
electron  levels  in  spherical  effective  potential. 

This  'work  was  supported  in  part  by  grants  from  the  Petroleum  Research 
Fund,  administered  by  the  ACS  and  the  Office  of  .Naval  Research  Orant 
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Figure  1 


FIGURE  CAPTION 


Electrical  polarizability  as  a  function  of  the  :luster  size.  N 
is  the  number  of  a  tens  in  the  cluster,  «  the  pel arizabi 1 i ty 
of  the  cluster  and  ji  the  polarizability  of  a  free  atom.  The 
black  dots  are  theoretical  results  of  the  jolliun  model  and  the 
open  circles  experimental  average  pol arizabi 1 i ties  for  ref.  1. 
The  filling  of  energy  levels  in  the  jell  inn  model  are  marked  by 
arrows.  Note  that  different  values  x.  are  used  for  scaling  the 
theoretical  (om=131  a;;)  and  experimental  (  ^ :  =  1 5 ? v. )  results. 


